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An attractive proposal for the hidden order (HO) in the heavy electron compound URu 2 Si 2 is 
an itinerant multipole order of high rank. It is due to the pairing of electrons and holes centered 
on zone center and boundary, respectively in states that have maximally different total angular 
momentum components. Due to the pairing with commensurate zone boundary ordering vector the 
translational symmetry is broken and a HO quasiparticle gap opens below the transition temperature 
Tho- Inelastic neutron scattering (INS) has demonstrated that for T < Tho the collective magnetic 
response is dominated by a sharp spin exciton resonance at the ordering vector Q that is reminiscent 
of spin exciton modes found inside the gap of unconventional superconductors and Kondo insulators. 

We use an effective two-orbital tight binding model incorporating the crystalline electric field effect 
to derive closed expressions for quasiparticle bands reconstructed by the multipolar pairing terms. 

We show that the magnetic response calculated within that model exhibits the salient features of 
the resonance found in INS. We also use the calculated dynamical susceptibility to explain the low 
temperature NMR relaxation rate. 

PACS numbers: 74.20.Rp, 75.30.Mb, 75.40.Gb 


I. INTRODUCTION 

The continuing investigation of 5f heavy fermion com¬ 
pound (7 ~ 170 mJ/mol K^) URu2Si2 is mostly moti¬ 
vated by its exotic low temperature phases. Firstly a 
much discussed ’hidden order’ (HO) state appears below 
Tho = 17.5 K and at considerably lower temperature 
Tc = 1.45 K unconventional chiral d-wave singlet super¬ 
conductivity (SC) appears. In the p-T phase diagram su¬ 
perconductivity is embedded within the HO phase. We 
shall not discuss here the long history of that subject but 
refer to the review in Ref. [TJ 

In the early theoretical work on thermodynamical 
properties the HO phase was interpreted in terms 
of localized and crystalline-electric-field (CEF) split 
5f stated and their localized multipolar degrees of 
freedoirP. However, it has subsequently b ecom e evident 
from comparison of ARPES experiments^ and local 
density electronic structure calculationS^ that the 
itinerant 5f electron band picture is a more appropriate 
starting point. The LDA calculations show that the 
heavy Fermi surface in URu2Si2 consists mainly of 
electron sheets around the F-point and hole sheets 
around the body centered tetragonal (bet) Z-point 
(0,0, ^) consisting mostly of j = 5/2 5f orbitals. Due 
to hybridization and Coulomb interactions the effective 
masses of the main Fermi surface (FS) sheets observed 
in de Haas-van Alphen (dHvA) and Shubnikov-de 
Haas (SdH) experiments are enhanced by a factor 
m* /rub ~ 7 — 8 as compared to the band masse^^. The 
sheets around F and Z are nested by the commensurate 
wave vector Q = (0,0,^) which is also the ordering 
vector of the HO phase. This translational symmetry 


breaking is suggested by ARPES features that indicate 
a downfolding of the Fermi surface from bet to a simple 
tetragonal (st) structure connected with gap opening 
in the nesting regiorP^. In addition HO breaks the 
tetragonal in-plane symmetry from C 4 to C 2 as observed 
in torqu^ and cyclotron resonanc^^ investigations. 
There is also direct evidence from x-ray diffracti on of 
a corresponding reduction in lattice symmetr^ bdi^ l. 
Furthemore analysis of NMRP^ and fxSI¥^ experiments 
suggest the breaking of time reversal symmety in the 
HO phase. Because the tetragonal symmetry breaking 
can only be achi eved by a doubly degenerate HO 
representatiorP^^R^ the i?_ representation is singled out 
as best candidate for HO. We stress, however that nu¬ 
merous alternative theoretical models of both localized 
or itinerant character have been proposed for the HO 
in URu2Si2. A partial list of these models and their 
symmetries has been compiled in Ref. |T3 (Supplement) 
and RefITCl 

In the itinerant picture the HO parameter may be con¬ 
structed as an electron-hole pair condensate of the nested 
Fermi surface states. Because the antiferroic (AF) -type 
Q connects two nested sheets with very different total 
angular momentum componentpRTI = m, M' the HO 
parameter is not of the usual spin- or charge density wave 
but has a high multipole rank r = max\M — M'\ due to 
their pairing. This has been explicitly shown in an ab- 
initio based multi-band interaction modeP^ where indeed 
the leading instability is of type with rank r=5 (’do- 
triacontapole’). Therefore the e-h pairing of states which 
contain jz = ±5/2 components is realized leading to the 
high-rank multipole character of the HO. Alternatively 
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to using the ab-initio band structure one may devise a 
simplified effective two-orbital model for the quasipar¬ 
ticle bandJi^ which can be used more easily to study 
the effect of HO and associated gap opening on vari¬ 
ous physical problems like quasiparticle interference^, 
thermodynamic?^, transport quantities and in particu¬ 
lar inelastic neutron spectroscopy (INS) which is investi¬ 
gated here. 

The latter has demonstrated that within the quasi¬ 
particle charge gap A//o = 4.1 meV that opens in 
the HO phase and is confirmed in STM experiment?^ 
a pr onoun ced and well-developed spin resonance 
evolve?32ll at the ordering wave vector Q = (0,0, f) 
(which is equivalent to (^,0,0) in the bet BZ). It 
has an order-parameter like temperature dependence 
of resonance frequency and intensity, where the former 
reaches ojr = 1-7 meV for T Tho corresponding 
to LOrlAno = 0.41. The spin resonance is a com¬ 
mon phenomenon in the gap of, e. g. heavy fermion 
(HF) unconventional superconductor?2Hlll^ but it may 
also exist in the hybridization or hidden order gap 
of non-superconducting Kondo compound?211231 (for a 
review see Hef lM|l . As in the case of Kondo insulator 
YbBi 2 the spin resonance in URu 2 Si 2 shows some 
considerable dispersion when moving away from the 
commensurate Q into the Brillouin zone. There is 
a further inelastic excitation at an incommensurate 
wave vector Qi = (1.4—,0,0) known already much 
earliei!2SEIl which, though broadened, exists even above 
Tho and therefore is not directly linked to HO param¬ 
eter and gap opening. It has been discussed earlie?^ 
and will not be addressed in the present HO context. 
We focus exclusively on the low energy spin resonance 
at commensurate Q which appears only below Tho and 
is clearly associated with HO gap formation because, 
unlike the excitation at Qi it vanishes for p > pc = lA 
GPa in the AF phas?^^. 


ture functions. In Sec. VH we briefly discuss interpreta¬ 
tion of NMR within this approach and Secs. |VI|VHI| give 
the discussion of numerical results and conclusion. 


II. TWO ORBITAL MODEL FOR HIDDEN 
ORDER AND HEAVY QUASIPARTICLE BANDS 

The 5f quasiparticle bands in URu 2 Si 2 are the result 
of hybridzation of narrow 5f bands exhibiting strong 
Coulomb interactions with wider conduction ba nds. 
Their enhanced effective masses (m* = 8 — 25too)P^ can¬ 
not adequately be described by band structure calcula¬ 
tions (mb < 8mo)P although the latter give a Fermi sur¬ 
face sheet topo logy in agreement with dHvA and SdH 
experiment?^. Conceptually mass enhancement can 
be obtained within the Anderson lattice model in con¬ 
strained mean field approximatioiP^ which also leads to a 
narrow indirect hybridization gap close to the Fermi level, 
equivalent to the Kondo temperature (Tr- ~ 11.1 m eV or 
129 K) that can be seen in STM experiment?SIlIl Here 
we are primarily interested in the very low energy mag¬ 
netic response of the heavy quasiparticle w ~ Aho =4.1 
meV which is considerably smaller than the Kondo tem¬ 
perature Th = 2.7 A HO- Therefore we do not start from 
an Anderson lattice type model with hybridizing narrow 
5f and wide conduction bands. We rather model directly 
the heavy quasiparticle bands and their Fermi surface by 
a simple ef fecti ve tight binding model. Ab-initio LDA 
calculation?^^ have shown that states at the Fermi sur¬ 
face consist mostly of jj-coupled total angular mometum 
5f states \j = |,M) with M = ±|,±| . Therefore the 
heavy bands, which have twofold Kramers- pseudo spin 
degeneracy may be directly described within an effective 
two-orbital tight binding model due to Rau and Ke?^ 
The basis states and pseudospins of this model are de¬ 
fined by 


The aim of our investigation is to understand and sim¬ 
ulate the resonance formation within the itinerant mul¬ 
tipolar A_(rank 5) HO scenario motivated above and 
based on the two-orbital model. We show that the prop¬ 
erties of the commensurate resonance found in INS can 
be described and understood semi-quantitatively within 
that itinerant model approach for URu 2 Si 2 by calculating 
the collective magnetic response in RPA approximation 
based on the four reconstructed quasiparticle bands in 
the HO phase. We also show that the hitherto not un¬ 
derstood NMR relaxation behavior in the HO phase may 
be qualitatively reproduced by the present model as the 
result of opening of a global spin gap. 


In Sec. |TT] we recapitulate the basic ingredients of the 
model and give the closed solution for the four HO quasi¬ 
particle bands. The total angular momentum operators 


in itinerant basis representation are given in Sec. HI 


Then we calculate the bare and collective RPA magnetic 
response functions in Secs. |IV| and |Vj respectivly which 
leads to the magnetic excitation spectrum and INS struc- 


( h± \ - f cos6» sin6> \ f /±f V 
\f 2 ±) V -sin^ cos9 J \ f^3 ) ’ 


( 1 ) 
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where A = x,y,z. Here a = 1,2 is the orbital, ± their 
Kramers pseudo spin degeneracy index and is the 
pseudospin in terms o f the orbital creation and annihila¬ 
tion operator j^^b7 | 40 | mixing angle 9 results from 
the CEF potential and is determined by the uniaxial 
anisotropy of the susceptibility (Sec. VI). In this basis 
the kinetic energy may be written as 


Hq — (Aik/j£,j(./lCTk + A2k/2o.k/2tTk) + 
kCT 

[7?k(/i+k/2-k — / 2 +k/l-k) + H.C.] , 
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FIG. 1. (Color online) (a) Fermi surface of the paramagnetic phase downfolded with Q to the st BZ of the HO phase consisting 
of two closed interpenetrating sheets, b) FS cuts in the = 0 plane for paramagnetic (PM) phase {(j> = 0; blue dashed line) 
and HO phase ((j> = 1 with red solid line). Momentum range in (b) is given by —rr < ki < n. (c) Fermi surface sheets of 
URu 2 Si 2 in the E-{1, 1) HO phase {(p = 1) in the st BZ. The FS is reconstructed in k-space regions connected by the nesting 
vector Q = (0,0,1) and breaks up into four smaller sheets without C 4 symmetry (for Ife^l > 0). 


where the model parameters of Refs. [TH] and [12] are im¬ 
plied which are also given in Appendix Following 
Refs. ITTUIOI the HO is a particle-hole condensate of mul¬ 
tipolar rank 5 which is breaking the following tetrago¬ 
nal space group symmetrie^i^ i) translational sym metry 
due to AF structure of HO with Q = (0,0, ii) ro¬ 

tational in-plane symmetry C 4 which is reduced to just 
C 2 symmetrjl^ iii) time reversal symmetrjff^. It is de¬ 
scribed by a doubly degenerate representation of the 
order parameter </> = {(px,(py). Using the pseudospin 
representatiorP^E^Eni (j = gj^e) rank-5 ’dotri- 

akontapole’ is defined bjff^EH 


= ^{3^12 + ^y- (0 = 

(4) 

In the HO phase one has to add the molecular field term 
to Ho {k = l/(2v^) = 0.35): 

■ Xl(/^k'^/ 2 k+Q + /Ik'^/ik+q) + H.c. (5) 

k 

The total mean field Hamiltonian, including the HO 
molecular fields may be written as 

H = Hq ~\~ H(j) = ^ ^ ^k — ^ak 0 ^fek- ( 6 ) 


It consists of two = a,b is the effective Kramers degen¬ 
eracy index) 4x4 blocks due to combined time reversal 
and translational symmetry T O tq . Defining the spinor 
basis 


V'lk = ifi 


1+k 

t 


ft ft 
> J 2 -k> J] 


l2-k’ 
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l-l-k+Q’ ■'2-k+Q 




— t ft 
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f)k ~ (/l-k’ /2+k> /l -k+Q> J2+\i+Q/ 

~ (l^ 6 kl’V^&k 2 ’V’Ls’^hk4)’ 


(7) 


and using /ik = (/ik — w/) = haia © h^k the blocks h-yk in 
corresponding order to Eq. ([^ are given by: 


^ak — 


hbk — 


/ Aik 

Dk 0 

A_ \ 


Dl 

A2k A+ 

0 


0 

A- Aik+q 

TAk+q 

5 

V 

n D* 

^k+Q 

A2k+q y 


( Aik 

0 

1 

A+ 

-Dk 

A2k A_ 

0 

0 

A+ Aik+q 7A^+q 


y A_ 0 —Dk+q A 2 k+q / 


( 8 ) 


Here A^k (a = 1,2 orbital index) and Dk are defined 
in Appendix Furthermore A± = with 

K = l/2-\/2 = 0.35. We note that hak t-t ^hk under 
transformation Dk —t D* = —D^ and A± —t Aj- = A^;. 
Although the A_ HO breaks time reversal, the product 
with translation tq (k —>■ k -|- Q) is still a preserved 
symmetry. Therefore the four (12 = 1 — 4) quasiparticle 
bands obtained from solving \h-yk — wl| =0 still have an 
additional twofold quasi-Kramers (7 = a, b) degeneracy. 
They are given by (12 = ±, 1, 2; = 12 = 1 — 4 

^ik — “ 

±^A^^+Tl + n^\cp\^ ± 2^A^^Tl + 

(9) 

where F^ = + |7IkP- On the l.h.s. ± corresponds 

to the second (under the square root) and 1,2 to the 
first (outside the square root) ± on the r.h.s., respec¬ 
tively. Expressed in a simplified way, Aj^ is the splitting 
of CEF orbitals 1,2; Ak is an intra-orbital and Dk an 
inter-orbital (~ ti 2 ) hopping energy, respectively (Ap- 
pendix[A| . Furthermore 4> = \cf)\ = {(1)1+(py)^ is the order 
parameter amplitude with cp = (p/\(p\ = {(j)x,(i)y)- Using 
Dk = D'^ + iD'^ we introduce the azimuthal function Ck 
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which leads to the breaking of fourfold (74 symmetry in 
the HO phase. It has the general form 

Ck = ^(|i?kp-^k) = 

i [D'^^ + - {D'J, - D'i^yf]. 

( 10 ) 

Due to Dk_|_Q = — Dk it also satisfies Ck+q = Ck- Here 
we consider only the phase with (px = 4>y and 

then it simplifies to 

Ck = =32^12 cos^ sin^ “fcj, sin^ 

( 11 ) 

which shows manifestly the fourfold symmetry breaking 
under (kx,ky) —>■ {ky,—kx) ■ However, the quasiparticle 
bands exhibit the translational symmetry £i/k+Q = Si/k 
of the HO phase. 

The unitary transformations Lt^k that diagonalize 
are defined by {{i', i/' = 1 — 4) HO quasiparticle band in¬ 
dex and 7 = a, 6 effective Kramers degeneracy of bands): 

^-yk — ^7k^7kbCyk’ — £z/k^z/i^S ( 12 ) 

Here Cj^k is degenerate in 7 = a, b. The columns of (/(Jk 
are the eigenvectors of the block matrices /i..yk- The cor¬ 
responding HO quasiparticle operators are defined by 

'07k ll7k'07ki 07k ^7k07k' (13) 

They satisfy canonical commutation relations 

{' 0 (- 7 k>' 0 (£ 7 'k'} = ^i'!^'^ 77 ' 4 k' and lead to a diago¬ 
nal representation of the Hamiltonian given by 

H = Y1 £!^k0’J)yk0('7k- (l-l) 


The Fermi surface of the disordered phase {(p = 0) 
corresponding to the quasiparticle bands in Eq. ^ is 
shown in Fig. downfolded to the st Brillouin zone 
of the HO phase. It consists of two interpenetrating 
electron and hole surfaces whose k^ = 0 cut is shown 
in Fig. by dashed lines. The points of crossing (in 
the downfolded or projected picture) of both sheets are 
connected by the nesting vector Q and therefore when 
HO appears {p > 0) these regions become gapped and 
the Fermi surface turns into dissected petal-like shapes 
shown as full lines in Fig. [^. This is also visible in 
the 3D FS presentation of the HO phase in Fig. [^3. It 
is obvious that large parts of the Fermi surface are re¬ 
moved in the HO phase. These model FS sheets cor¬ 
respond well to the (3— band sheets from LDAI^ calcu¬ 
lations and dHvA and SdH experiment ^311 as well as 
ARPEpSl results. They are also the heaviest sheets 
with m*p ~ 25mo. From the shrinking of the FS in 
HO phase one expects that a pronounced suppression 
of the DOS at the Fermi level (w = 0) appears. This is 
shown in Fig. in comparison to the disorderd phase. 



FIG. 2. (Color online) (a) Density of states (DOS) 
showing the evolution of HO gap for small ui. (b) Tem¬ 
perature dependence of DOS at the Fermi energy {lj = 
0). Here we used a mean field-type interpolation p{T) = 
<)> tanh( 1.76 -^/Tho/T ~ 1) for the temperature dependence of 
HO parameter. 


Then, using a mean-field type interpolation for the tem¬ 
perature dependence of the order parameter according to 
p{T) = (^tanh(1.76-Y/Tf/o/T — 1) the DOS at the Fermi 
level (w = 0 ) is progressively reduced to about one half 
of the original value (for p = 1) when temperature drops 
below Tho (Fig.[^). 

The breaking up of the large FS sheets into smaller 
petals and their downfolding to the F point in the HO 
phase was found in ARPES experiments^. The reduc¬ 
tion of the DOS explains the observed entropy loss and 
suppressed 7 -value in the HO phasS^. 


III. THE MAGNETIC MOMENT OPERATORS 
IN THE TWO-ORBITAL MODEL 


The dynamical susceptibility which appears in the INS 
structure factor and cross section contains the total an¬ 
gular mometum operator of 5f states. Therefore we need 
its representation in terms of single particle creation and 
annihilation operators of j = | eigenstates. It is defined 
by 


Jx=J 2 


(15) 


MM' 


with the explicit matrix elements given in Appendix 
This leads to 

=^x{flh+flh)^ + Kfl.f-|+fKf J., 

2 ^ 2 ^ 2 ^ 22 

Jl = - h - ft h )i + *m(/1 5 /- 3 - /! 3 / 


=|(/i/i - /!./-i)* + |(/l/i - /L/-f).- 

Z 2 ^ 2 ^ Z 2 2 ^ 


( 16 ) 
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T(K) 


FIG. 3. (Color online) Temperature evolution of the different 
tensor components of static susceptibility Xo{q = 0,tu = 0): 
(a,b) XX and 22 : components exhibit suppression below Tho- 
Note the scale difference due to larger magnetic moment ma¬ 
trix elements for z-direction (c) xy component vanishes above 
Tho because of tetragonal symmetry and becomes finite be¬ 
low due to fourfold symmetry breaking induced by E-{1, 1) 
HO phase. 


The magnetic moment operators can also be written in 
pseudo-spin components at site i {a, a' = 1,2 or¬ 
bital index, A = x,y,z cartesian index) or in terms 
of their Fourier components. They are obtained from 
Eq. as 


= — V 

/ V j (y.<T\L^ era' J ( 


Q;fTk^o-o-'/Q!'o-'k+q- 


(17) 


Using the transformation of Eq. ([^, inserting into 
Eqs. (16) and performing the Fourier transform we gelP^ 


=v^[cos 26{S^^ + + sin 26{S'^^ - S ^^)], 

jq =75 [cos 26{S^y + 52 T) + sin 2e{S'^f - S^D], 
Jf =4[cos29iS'^^ - - sm20{S^^ + S^^)] 

+ {S^1+S^^). 


The total angular momentum components in Eq. (18) 
may also be expressed as bilinear forms of spi nor opera¬ 
tors -tpl = of Eq. ([^ by using Eq. (17). Then 

we have to restrict summation over k to the nrst st BZ 
of the HO phase and obtain 


~ ^ V'kEA^’k+q 

(19) 

Here the signs sx are given by s^, = 1 and Sy = Sj = — 1. 
The F^ F^ matrices are given by 


F° =ToKxfJ. cos 26-, F^ = ToKzTsm29, 

F° =To{iKy){iy) COS29-, FJ, = ToKo(iAi) sin20, 

r° =ro(2cos20Ko + = —To{iKy)2sm29. 


( 20 ) 


where Kx,Ky,Kz are the Pauli matrices and kq the unit 
in orbital (a = 1,2) space while tq is the unit in HO 
Nambu space. The explicit r)[’ matrices are presented 
in Appendix [B| 


Finally we have to transform to the spinors 
(7 = a, b) which diagonalize the Hamiltonian in the HO 
phase. This is achieved by the unitary transformatiorfl^ 
of Eq. ( [I^ resulting in 

= (21) 

The transformed moment matrix elements which are now 
momentum dependent are given by (A = x, y, z) 


( tlakPiC/lk+q \ 

^ SAC/.kPit/ik+q SAt7,kr°,Hlk+q J 

aa Wab 

^ Akq ^ Akq 
T^'ba 'p'bb 
^ Akq ^ Akq 




( 22 ) 


Here again Sx = l,Sy = Sz = —1. The essential matrix 
elements for calculating the bare magnetic response func¬ 
tions y7'(Hj*<A’n) in the HO phase are then contained in 
the matrices F^[[q (l,l' = a, b) where F°’^ are defined in 
Eq. (B3) and is the unitary matrix that diagonalizes 
the 4x4 block Hamiltonian matrices h-yk of Eq. (1^. 


IV. THE BARE MAGNETIC SUSCEPTIBILITY 
IN THE HO PHASE 


The physical magnetic moment is then obtained through 
H = ( 5 Mb)J with g = 


The bare dynamic magnetic susceptibility (A, A' = 
a:, ?/, z are cartesian matrix indices) is obtained from the 
Fourier transform of magnetic moment Green’s function 
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FIG. 4. (Color online) Bare (noninteracting) dynamical sus¬ 
ceptibility Xo^(Q)<^ at the bet Z-point (1,0,0) (r.l.u.) below 
(0=1) and and above (0 = 0) HO temperature: (a) is the 
imaginary part, and (b) shows the real part. The singular be¬ 
haviour in the low energy region is due to the HO gap opening 
and leads to the resonance appearance in Fig. 


Changing the summation in Eq. (25) from k into —k and 
using inversion symmetry we can finally write 




l/'k-f-q 


- /l/k 


(27) 

Using the explicit component representation of momen¬ 
tum operator of Eq. (1^ according to 


= ^ U (28) 


in Eq. (26 1 we then obtain the matrix elements 


^AA'(kq)= (l^:k + q + K))p 

77 'K 

_ J_ ('-p'77''ikq+K*/p'77'\kq+K 

/ j A )vv' X' Juv' 

77 'K 


(29) 


according to 

xo^'(q,w) = 

1 (nk| J^|n'k')*(nk| J^,|n'k') 

'n ^ 


k,n,k'n' 


^nk ^n'k' ^ 


(/n''k' /nk)- 


Here n = {v, 7) is a double index with u = 1 — A denoting 
one of the four HO bands and 7 = a, & is the Kramers 
degeneracy. One can sum over the latter only in the 
matrix elements since £„k = ei/k- To obtain the correct 
singular behaviour of (q, w) close to the HO gap it 
is essential to include the reconstruction of Bloch states 


by the HO parameter in the matrix elements of Eq. (23). 


They are nonvanishing when the sum of wave vectors 
corresponds to a reciprocal lattice vector K of the simple 
tetragonal (st) BZ of the HO phase: 

(nk| J^jn'k') = (n'k'| J^'*|nk)* 

= ^J_k+qs-k'+K.om™(kk'). (24) 

K 

Inserting this into Eq. ( [2^ we obtain 

xo^'(q>^) = 4 E E 


77'K 


Kir''(k,k-q-K; 


|2 /y'k—q /jyk 


(25) 


where we used the Kramers degeneracy e„k = £i/k and 
periodicity of bands in the HO phase £i^k-i-K = ^i/k- Ex¬ 
plicitly, we have 


^ Kir'"'(k,k-q-K))|2 = 


77 'K 


X! ('^Tkl J^lz^Vk - q - K)*(z/ 7 k| J^,|j^' 7 'k - q - K). 

77 'K 

( 26 ) 


For the HO phase with wave vector Q = (^,0,0) only 
K = 0, Q are independent because the matrices are 
periodic with 2Q and the energies with Q. Then the sum 
over K contains only two terms leading the final explicit 
matrix elements: 


(23) Mj:!(;(kq) = E[(rl""')L^'*(r;'("')::^' 


77 


(30) 


^(j,y)kq+Q*(rU7')kq+Q]^ 

in terms of the transformed moment matrices of Eq. 

They fulfil the periodicity (kq-f Q) = (kq) of 
the HO phase st lattice. Together with the periodicity of 
the energy bands this means Xo'^ (q+Q, w) = Xo^ (qj^;). 
For the diagonal cartesian susceptibility components the 
above matrix elements simplify to 

MC!('(kq) = 

[l(r;")::.^.p + + (q-> q + Q)], 


E 


with the explicit T^ matrix elements 


l(rE)::^f = l(c/ 7 kr^At 74 +q). 
l(rE)L^f = I(c/7kriu4+q). 


(32) 


and similar for wave vector q+Q. Here we used the con¬ 
vention 7 = 6,0 for 7 = Q, 6 respectively in the nondiago¬ 
nal parts. The Eqs.( 27|3'2 |, together with HO bands e^k, 
the HO transformation matrix C^k and the T^’^ rnatrices 
of Eq. (B3) are the ingredients to calculate the diagonal 


susceptibility elements. The no ndiagonal (xy-type) ones 
of general case in Eqs. ( 27|30 ) can only appear in the 
(trigonal) E_(l,l) HO phase as induced elements and 
therefore will be quite small. 
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FIG. 5. (Color online) (a) Collective magnetic {zz, _L) ex¬ 
citation spectrnm: the imaginary part of the collective RPA 
susceptibility at the bet Z point versus energy. The position 
of the pronounced peak determines the (almost T- indepen¬ 
dent) resonance energy in the HO phase. Here and in 
following hgures J_i(Q) = 0.156to = 1.04 meV (b) Evolu¬ 
tion of imaginary part of the RPA susceptibility at resonance 
energy uir{T) in q-scan around bet Z- point for different tem¬ 
peratures. 


V. THE HO-RPA SUSCEPTIBILITY AND 
APPEARANCE OF COLLECTIVE SPIN 
EXCITON RESONANCE AT Q = (1,0,0) 

The bare magnetic response will be enhanced by effec¬ 
tive exchange interactions between quasiparticles in the 
HO phase. This is also suggested by the proximity of the 
AF phase that exists above a small critical pressure. If 
the formation of heavy quasiparticles is described by the 
constrained (slave boson) mean-field approximation of an 
Anderson-lattice type model these quasiparticle interac¬ 
tions are caused by ffuctations beyond the mean field 
leveP^. However their momentum structure is too singu¬ 
lar and therefore we use here a phenomenological form of 
the exchange interactions described by 

= JAA'(q)J^V. (33) 

qAA' I 




1 ({^-J\\xl^)xr 

i:i||(q,w) (y (i-^iixDxr 


Here the determinant is obtained from 

D|,(q, cu) = 1-J|| (q)(xS"+Af(")+J|| (q)^(xr Xr-Xo^xD- 

(38) 

For vanishing non-diagonal {xy, yx) susceptibility ele¬ 
ments x']jp^(q, w) is proportional to the unit matrix and 

Xppyi(q, w) is obtained from Eq. (361 by replacing T—>-11. 
The dynamical structure function which is proportional 
to the INS scattering cross sectioil2iH2l of URu 2 Si 2 is 


Generally, within the RPA approach the effective ex¬ 
change function J(q) must have its maximum at the wave 
vector where the main resonance peak appears, which is 
the zone boundary (ordering) vector Q = (1, 0,0) (in this 
section r.l.u. units ^ of INS are used). The q depen¬ 
dence around Q may then be used for fitting to the range 
of the resonance dispersion. In URu 2 Si 2 the dispersion is 
well localized at Q with only weak and steep dispersive 
features close to itP2 as discussed below. The numerical 
results will show that this is dominated by the q- depen¬ 
dence of the bare susceptibility. Therefore we can use the 
simplest approximation of a q-independent J(q) = ^(Q) 
for the exchange. The numerical value of the interaction 
parameter J(Q) is then obtained by requiring the peak 
position in the RPA spectrum to be close to the exper¬ 
imental resonance position. In the tetragonal structure 
exchange function J(q) and susceptibility ;yo(qjw) are 
uniaxial tensors J(q) given by 

/J,l(q) 0 0 \ 

J(q)= 0 J,|(q) 0 ; 

V 0 0 J^(q )j 

/xr(q,^) xS"(q>^) 0 

xo(q,w)= 0 

V 0 0 A^(q,w) 


(34) 


The non-diagonal elements in the susceptibility matrix 
may appear for the £'-(1,1) HO phase but not for the 
£_(1,0), £-(0,1) and disordered phases. The collective 
RPA susceptibility then has the form 


x(q,w) = [i - xo(q,w)J(q)] ^xo(q,w) 

= ( xIpACq.^^) 0 
V 0 xkpAin.^) 


(35) 


where inversion leads to the final result 


Xrpa{^t^) 


1 - '/-L(q)x^(q,w) 


-1 

xo (q>^^)> 


(36) 


+ J\\x7xr (i-^iixr)xS" + ^iixW\ 
+ jyxrxr (1 - ^iixDxr+^iixo^xr 


(37) 


given by (/3 = (fcsT) ^): 

5'(q,w) = - 9AgA')^"T-XAA'(q,w). 

AA' 

(39) 

For small temperatures (/3a; 3> 1) the Bose-factor 
(1 — —>■ 1 and for large temperatures (/3a; 1) 

(l_e-/5“)-i ^ {(3uj)-\ 
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RPA expressions of Eqs. (36|37 ) we obtain: 




— Im 

TT 


Xo 

1 - Jj-Xo 


+ 


xr-Mxrxi^-x7xr) 


1 - Jwixo"' + xD + -^11 (x^x^ - Xo^xD^q’^' 


(41) 

We will later plot the two contributions separately for 
clarity. Because of the large uniaxial anisotropy of ma¬ 
trix elements (Appendix [b| the sum will be determined 
by the A {zz) contribution; for isotropic Jy = Jl only 
this channel develops the resonance peak. It is also 
useful to consider two special cases that may be realized, 
depending on the type of HO, i.e., whether in-plane 
tetragonal symmetry is preserved or broken and in the 
latter case whether (1,1) or (1,0)- type orientation of 
the degenerate order parameter is realized: 


i) Fourfold symmetry breaking but vanishing non di¬ 
agonal elements: 

This case corresponds to E_{10) or A_(01) with Xxx ^ 
Xyy and Xxy = Xyx = 0. Then we have 


5'(q,a;) = —Im 

TT 


Xo 


X^ 


L 1 -J_lXo 


(42) 


FIG. 6. (Color online) (a) Contour plot of the imaginary part 
of the bare susceptibility Xo^ A to tetragonal plane, (b) the 
RPA susceptibility Xrpa (c) corresponding plots of the imag¬ 
inary part of the bare and (d) RPA spectrum of {yy, ||) dy¬ 
namical susceptibility with a large J||(Q) = O.Qlto = 5.8meV. 
For istotropic J||(Q) = J±{Q) no resonance appears in this 
channel. Here we set T = 3K, in the HO phase with (j) = 1. 
The resulting resonance peak in (b),(d) is at uJr = 0.18 meV. 
The dispersive V-shaped tails are remnants of the quasi-lD 
features of particle-hole continuum due to the nesting. 


In INS investigations the scattering vectors for the 
bet structure are conventionally indexed by those of 
the simple tetragonal structure^^, i.e. by the carte¬ 
sian compo nents . Mainly the [100] direction has been 
investigatecPiEIl gofar where q = qq^,. Note that the in¬ 
plane zone boundary vectors Q = (1,0,0) and (0,1,0) 
are equivalent to the Z-point vector Q = (0,0,1) of the 
bet structure which is also the ordering vector of the 
antiferro- type HO. This is due to the fact that adjacent 
bet BZ’s are shifted along [001] directiorP^. For the [100] 
momentum direction we then obtain (i.f. /3a; ^ 1 or 
0 ): 


If the fourfold symmetry breaking is absent Xo^ = Xo^ = 
Xq. Then the above expression is formally the same as 
for the disordered phase {(p = 0) above Tho- 
a) Fourfold symmetry breaking through finite and equal 
non-diagonal elements: 

This case corresponds to i3_(ll) or E-{ll) with Xxx = 
Xyy = Xo and Xxy = Xyx = Xo leading to 


^(q, uj) = —Im 

TT 


Xo 


Ll - J_LXo 


A 


xl - J\\{xf - xl) - 

(1 - -^llXo)^ - J^iXl' 


q,u; 


(43) 

If the non-diagonal element xo induced by i3_(ll) HO is 
negligible this reduces again to the simple RPA expres¬ 
sion 


iS'(q, uj) = —Im 

TT 


Xo 


Xd 


Ll-J^Xo 1-J\\xl 


(44) 


which is formally identical to the disordered phase but 
Xo(q, a;) now contain the effect of HO band reconstruc¬ 
tion. For nonzero temperatures the above expressions 
have to be multiplied by (1 — 


5'(q,a;) = -/TO XRPAi^l^^) + XkpAi^^^^) 

TT L 


(40) 


VI. DISCUSSION OF NUMERICAL RESULTS 
FOR THE MAGNETIC EXCITATION 
SPECTRUM 


which means that only moment fluctuations A to the Now we discuss the numerical results for the spin 
[1,0,0] direction contribute to ^(q, w). Using the explicit dynamics in the two-orbital model of URu 2 Si 2 . First 






























9 




FIG. 7. (Color online) 3D plot 
of the imaginary part of the per¬ 
pendicular component of (a) bare 
(noninteracting ) Xo^ s-nd (b) 
RPA susceptibility Xrpa in hid¬ 
den order phase with with 0=1 
and T = 3K (same as Fig[^). In 
the latter main resonance peak 
and V-shaped dispersive features 
at higher energies can be identi¬ 
fied. 


we focus on the bare noninteracting susceptiblility. The 
static, homogeneous (q = 0 ) tensor components are 
shown in Fig. as function of temperature. The large 
uniaxial anisotropy Xzz/Xxx — 6.1 which is due to the 
CEF is comparable to the experimental on^^. It is 
used to determine the CEF mixing angle 9 ~ O.SIStt 
in Eq. ([^. Below Tho Xxx = Xyy as well as x^z 
show a considerable depression caused by the HO gap 
opening. Due to different selection rules for and 
Jz- operators for band states their behaviour around 
Tho is distinct. On the other hand the non-diagonal 
in-plane susceptibility Xxy which has to vanish above 
Tho becomes finite in the HO phase. This is the reason 
for the a ppea rance of twofold torque-oscillations in the 
HO-phase^^. 


mately reproduced. The latter is at ~ O.lSto = 1-2 
meV (to = 6.7 meV) or u}r{T = 0)/A//o — 0.29. This is 
in reasonable agreement with the experimental results 
Wr{T = 0)/Aijo — 0.41 where Wr = 1.7 meV^ and 
/S.HO =4.1 me\^. There is almost no temperature 
dependence of the resonance position once it exists. This 
is also in agreement with experimenlPlI^. On the other 
hand no clear indication of spin gap behaviour at the 
lowest u) precisely at the resonance vector Q is obtained. 
There is also a smaller side hump at larger energy due 
to the higher energy singularity in xo(Q,w) (Fig. |^. 
Alternatively we show a q-scan along {q, 0, 0) direction 
for constant uj = uJr{T) at various temperatures. Again 
the resonance peak clearly grows below Tho- 


The bare dynamical susceptibility (the Lindhard func¬ 
tion) at the Z-point is shown in Fig. The HO gap 
opening produces singularities in the response around 
the gap threshold oj ~ ^ho = </'/'\/2- These singularities 
are responsible for the resonance appearance in the 
collective response function. To obtain them one must 
use the reconstructed eigenstates and matrix elements of 
Eqs. ( 22|30 ). Note there is an additional peak behaviour 
at higher energies uj ~ 2 to which is connected with the 
van Hove singularities of the band structure itself and 
therefore persists above Tho when 0 = 0 . 


The imaginary part of the collective RPA susceptibil¬ 
ity in the T (zz) channel, i.e., the magnetic excitation 
spectrum of the interacting itinerant 5f moments is 
shown in Fig. [^for the experimental [100] direction. The 
frequency dependence in Fig. at the bet Z point (100) 
shows the evolution of the spin resonance excitation out 
of the normal state spin fluctuation continuum (T=18 
K) when the temperature drops below Tho = 17.5 K. 
The interaction parameters (see caption) in the model 
are chosen such that the enhancement of peak intensity 
I{T = 0)/I{Tho) — 4 relative to the disordered state 
and the position of the resonance at ujr are approxi¬ 


These results may also be demonstrated in contour 
plots of the magnetic spectrum in the q,uj plane (Fig.|^. 
In (a) the spectrum of the bare dynamical susceptibility 
Xo^ is shown which exhibits the HO gap at w/to > 1 
and a V- shaped structure emerging from the bet Z- 
point {q = ^). Turning on the quasiparticle interaction 
leads to the spectrum of the collective RPA susceptibil¬ 
ity Xrpa shown in (b). Most of the magnetic intensity is 
now concentrated at the bet Z point resonance energy ojr 
(c.f. Fig .[^). But still an indication of the V-shaped dis¬ 
persion is visible. Qualitatively it agrees with the experi¬ 
mentally observed dispersion^^ around Z. One may ques¬ 
tion about its origin. In this model the HO is driven by 
nesting between the two FS sheets which contains states 
with high angular momentum component. The nesting 
means there is a quasi-lD contribution to the spectrum 
of the bare susceptbility. In such quasi ID-situation for 
low energies and small momenta with respect to nesting 
vector it consists of a very narrow particle-hole contin¬ 
uum with a dispersion v^pq' where q' = q — and v^p is 
the Fermi velocity in x-direction. The V-shaped ridges in 
(a,b) are their image. The spectrum of bare and RPA in¬ 
plane susceptibility {yy, ||) shown in (c,d) respectively de¬ 
pict a qualitatively similar behaviour as described above, 
although the scale is much smaller due to the uniaxial 
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matrix element anisotropy origninating in the CEF ef¬ 
fect. Therefore its contribution to S'(q, w) in Eqs. ( |40pT| ) 
is less important. Finally, the bare and RPA spectrum 
of magnetic excitations (zz, _L channnel) are presented in 
Fig.[7fi ,b , respectively in a 3D perspective plot where the 
described features of resonance formation and attached 
ridge-like dispersion can be seen even more clearly. 


VII. SPIN GAP FORMATION AND 
INFLUENCE ON NMR RELAXATION 


It is known from many examples that there is a dual re¬ 
lationship between spin resonance formation around the 
gap threshold (in this case HO gap) and a spin gap forma¬ 
tion at low energie^^ which are both seen in INS exper¬ 
iments. The spin gap formation also directly influences 
the NMR relaxation rate which is determined by low en¬ 
ergy spin fluctuations. If the latter open a spin gap the 
relaxation rate should strongl y decr ease. This was indeed 
found in the NMR experiment^^^l^in URu 2 Si 2 below the 
hidden order transition Tho- The normalized NMR re¬ 
laxation rate at NMR resonance frequency ujq is given via 
the dynamical RPA susceptibility by the relations 




/mx||(q,a;o) 

Wo 


,2 1 r-f™X||(q,Wo) 

9 L 


Wo 


Imx±(q, wq) 

Wq 


(45) 


where to is the energy scale of HO bands used previously. 
The relaxation rates are normalized to the Korringa rate 
(l/ri)o = 2xlAlj:kBT/{xeh)Hl where Ahf is assumed 
as a q- and axis- (||,-L) independent hyperfine coupling. 
The 7 „, 7 e are nuclear and electronic gyromagnetic ra¬ 
tios. To comply with previous convention (||,-L) denotes 
directions parallel and perpendicular to the tetragonal 
ab-plane (note this is opposite to conventions in Ref. \AS\ 
and fi^l) . The calculation of relaxation rate requires the q- 
dependence of dynamical susceptibility at resonance fre¬ 
quency Wq in the whole st HO Brillouin zone, not just 
along the Q-direction. The results are shown in Fig. 
for a frequency wq ^ he., much smaller than the 

HO gap. Below Tho a pronounced reduction is observed 
which shows that an overall low energy spin-gap is de¬ 
veloping in the HO phase although it is not localized in 
momentum space around the resonance vector Q because 
the HO gapping is rather incomplete. Qualitatively it 
is very similar to experimental results for ^®Si-NMR in 
Ref. 1151 We only present the (1 /Ti)_l rate. The corre¬ 
sponding (1/Ti)|| rate shows almost identical behaviour 
except for an enhancement factor due to the suscepti¬ 
bility anisotropy which is close to the experimental en¬ 
hancement (1 /Ti)||/(1/Ti)_l of about 4.d^. 



FIG. 8. (Color online) Temperature dependency of the 
normalised NMR relaxation rate, {Ti)±, at NMR resonance 
frequency wo ^ utr- Relaxation rate (ri)|| has similar T- 
dependence but is larger due to uniaxial anisotropy of mag¬ 
netic matrix elements. 


VIII. SUMMARY AND CONCLUSION 

In this work we have given a theoretical analysis of 
the spin-resonance phenomenon in HO phase of URu 2 Si 2 . 
The physical properties of this collective excitatio n has 
been well investigated before by INS experimentd^^Hl. 
It appears inside the HO charge gap Aho — 4.1 meV 
at an energy w^ = 1.7 meV at the commensurate wave 
vector Q = (^,0,0) of the bet Z-point and it exhibits 
a V-shaped upward dispersion. Such spin-resonance ex¬ 
citations are frequently observed inside the quasiparti¬ 
cle gap of unconventional superconductors like cuprates, 
pnictides and heavy fermion metals. More rarely they are 
found within the hybridization gap of Kondo insulators 
or the gap of hidden order compounds as in the present 
case. 

For a semi-quantitative explanation of this interest¬ 
ing many b ody e ffect in URu 2 Si 2 we used a previously 
investigate d^^F^ I two-orbital model of 5f electrons result¬ 
ing in two FS sheets with states that contain high an¬ 
gular momentum components. The hidden order is then 
taken as the rank-5 doubly degenerate E_ representa¬ 
tion which agrees with all observed broken symmetries. 
Within a mean held HO theory the quasiparticle states in 
the HO phase are reconstructed and a HO gap is opening 
leading to a vanishing of large parts of the Fermi surface 
and associated DOS suppression of heavy charge carri¬ 
ers. The reconstructed states lead to singular behaviour 
of the bare magnetic response function. This enables the 
RPA response function of interacting quasiparticles to 
develop a pronounced resonance peak below the HO gap 
at the nesting vector Q with a temperature evolution 
below Tho that is similar to the experimental one. In 
addition V-shaped dispersive features appear at higher 
energies which have also been observed in INS. Although 
there is no pronounced spin gap evolution at Q itself the 
momentum integrated low energy response seen in NMR 
develops such a suppression which is seen in the decrease 
of the NMR relaxation rate below Tho- 
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Embedded in the HO phase an unconventional SC gap 
of presumably chiral d-wave nature appear^^^SH. Its ef¬ 
fect on QPI was studied previousljff^. Because T^. ^ Tho 
by an order of magnitude the same is true for the super¬ 
conducting gap with respect to HO gap. Therefore we 
were not able to see a signihcant change of the magnetic 
RPA response function when including the superconduct¬ 
ing gap. This is also in agreement with INS finding^^. 
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Appendix A: Kinetic energy coefficients and para 
phase band structure 


Eor completeness and convenience we recapitulate 
the effective 5f- two band model for URu 2 Si 2 that is 
adopted from Ref. HH] and also used previously in QPI 
calculationJi^. The kinetic terms in Eq. ([^ are defined 
by the intra-orbital energies (a = 1,2 is the orbital in¬ 
dex): 

Aak =Kk + ^ak + ^sign{a)Ai2, 

=2i^(cos akx + cos aky) + cos ak^ cos aky — eo, 

.. ^ CL a, c, 

-^ak —COS COS COS -^kzf 

(Al) 

where we defined sign{a) = (—1)“ ^ and the inter¬ 
orbital hopping energy 

=ti2 [sin^(fca; -I- ky) - isin - ky)] sin|fcj;, 

(A2) 

To reproduce a realistic Fermi surface model with nest¬ 
ing electron- and hole- like pockets around the T and 
Z points of the bet Brillouin zone we use the following 
parameterJi^ The orbital energy splitting is A 12 = 3.5 
or A = O. 5 A 12 = 1.75. The nearest neighbor hoppings 
are ti = t 2 = t = —0.3, meaning orbital-independent 
= A^. Hopping elements to next and second nearest 
neighbors are = —0.87, t 2 = 0.0, f" = 0.375, f 2 = 0.25, 


/ 0 g cos 29 0 \ 

I /r cos 20 0 0 I 

0 0 0 ^cos20 ’ 

Vo 0 g cos 20 0 / 


respectively and the average orbital energy is —eo = 0.5. 
The inter-orbital hopping is |ti 2 | = 0.7. All energies are 
given in units of to = 6.7 meV corresponding to a total 
band widtll ^ Won — 12to = 80 meV obtained from tun¬ 
neling resultd^oisg] pgj. computation of quasiparticle 
bands in the HO phase it is useful to introduce (anti-) 
symmetrized quantities: 

Ak =2 + A^k) 

=2f'(cos akx + cos aky) + At" cos ak^ cos aky — eg, 

Ak =A + A4) 

=A -I- 25' [cos akx + cos aky) + A5" cos akx cos aky. 

(A3) 

Here t' = ^[t[ +^ 2 ), t" = |(f" -I- 1'^) and <5' = |(f( — ^ 3 ), 
5" = |(f" — The auxiliary functions above have the 
following symmetry under translation by the ordering 
vector Q- A^k+Q — ^ak implying also Aj|:_^q = Aj^ and 
Ak+Q = Aj^. On the other hand = — A^j^. and 

Dk+Q = —71k- 


Appendix B: Total angular momentum and 
susceptibility matrix elements for special cases 


Here we give the explicit form of the magnetic moment 
matrices in the basis of free single ion \j, M) states in the 
|j, ±|), |j, ±f) subspace that are needed to construct 


the total angular momentum operators in Eq. (161. Jx 


(A = a;, y, z) can be written as (index M, M' in the order 
f, f,-|, |), defining g = = 1.12: 


Jx - 


^ 0 /r 0 

g 0 0 0 

0 0 0 g 

\0 0 g 0 


X = 



\ 0 0 




0 

0 

-ifj, 

0 


(Bl) 

The susceptibility matrix elements for the normal 

state (</) = 0) are defined in Eq. (29) via the uni¬ 


tary transforms of the T^^ matrices in the space of 


four-dimensional '!/' 7 ,k, (7 = o-,b) spinors in Eq. (22). 
Explicitly, in each of ij^ak or 'i/j^k subspace they are given 
by (0 = CEE mixing angle): 


g sin 20 0 0 0 

0 —/i sin 20 0 0 

0 0 /i sin 20 0 

0 0 0 —/r sin 20 


(B2) 











12 


/ 0 

ifj, cos 20 

0 



/ i/isin 20 

0 

0 


1 —ifj, cos 29 

0 

0 

0 1 

1 ■ - 1 

' 0 

ifi sin 20 

0 

0 

0 

0 

0 

cos 20 j 


> 0 

0 

ifj, sin 29 

0 

V 0 

0 

—ifi cos 20 

0 ) 

' 1 

^ 0 

0 

0 

ifj sin 29 / 


(B3) 


r 


0 

Z 


/2cos20+i 0 

0 2 cos 20 

0 0 

\ 0 0 


0 

0 

2 cos 29 

0 


1 

2 

2 cos 


0 

0 

2 ^- 5 / 


0 —2 sin 20 0 

2 sin 20 0 0 

0 0 0 

0 0 2 sin 20 



(B4) 


Now we derive the susceptibility matrix elements 
for the special case I?k = 0 and without HO {(j> — 0). 
Then Hamilton block matrices of Eq. Q are already di¬ 
agonal and therefore trivially Ua,b = 1- The are 


therefore momentum independent and equal to the bare 
Fa matrices. Then the diagonal susceptibiltiy matrix el¬ 
ements of Eq. (32) reduce to the momentum independent 
expressions given by 


Mj:i' = 4[|(r°w|^ + |(riwp]. (B5) 


Written in concise 4x4 matrix form M\\ in {vv') indices 
we get for the three cartesian directions: 

Mxx = Myy = 4/x^To(sin^ 29kq -\- cos^ 29Kx), 

Mzz = 4to [(4cos^ 29 + + 2 cos 29Kz + 4sin^ 20^2,], 

(B6) 

where are the units in Nambu and orbital space, 

respectively and Kx,Kz are orbital Pauli matrices. In 
explicit matrix representation in -spaces (with order 
Ik, 2k, Ik -f Q, 2k -t- Q 51 ? 2 = orbital index) we have 


Mxx = M, 


yy 



0 0 \ 

0 0 

sin^ 20 cos^ 20 ’ 

cos^ 20 sin^ 20 / 


/ (2 cos 20-I-4sin^20 


M,, = 4 


0 

4sin^20 (2 cos 20—0 


0 

0 

dsin^ 20 


(2 cos 20 -I- 

4sin2 20 (2cos20-i)2y 


(B7) 


(B8) 


The sum of all susceptibility matrix elements fi\ = 
Si/i/' ^x\ is a measure of the effective moment in direc¬ 
tion A. They are independent of the CEE mixing angle 0 
because the latter causes just a rotation in the j=5/2 sub¬ 
space. We obtain /i^ = 64 and ^ix = ^J■y = = 10.03 

leading to an anisotropy ^z/^J■x=G■4^ independent of 0. 
This is approximately the anisotropy ratio of static sus¬ 


ceptibility. The blocks with zeroes in M\\ appear be¬ 
cause for (}() = 0 no mixing of states with momenta k and 
k-fQ is present. Finally we note that for ^ = 0 we also 
have Mx\' = 0 (A ^ A') because of preserved tetrago¬ 
nal symmetry. Therefore the y tensor in the disordered 
phase is uniaxial with Xxx = Xyy 7^ xL- 


* alireza@apctp.org 

^ J. A. Mydosh and P. M. Oppeneer, Rev. Mod. Phys. 83, 
1301 (2011). 

^ P. Santini and G. Amoretti, Phys. Rev. Lett. 73, 1027 
(1994). 

® Y. Kuramoto, H. Kusunose, and A. Kiss, J. Phys. Soc. 
Jpn. 78, 072001 (2009). 


^ J.-Q. Meng, P. M. Oppeneer, J. A. Mydosh, P. S. Rise- 
borough, K. Gofryk, J. J. Joyce, E. D. Bauer, Y. Li, and 
T. Durakiewicz, Phys. Rev. Lett. Ill, 127002 (2013). 

® C. Bareille, F. Boariu, H. Schwab, P. Lejay, F. Reinert, 
and A. Santander-Syro, Nature Communications 5 (2014). 

® P. M. Oppeneer, J. Rusz, S. Elgazzar, M.-T. Suzuki, T. Du¬ 
rakiewicz, and J. A. Mydosh, Phys. Rev. B 82, 205103 
( 2010 ). 













13 


^ H. Ohkuni, Y. Inada, Y. Tokiwa, K. Sakurai, R. Settai, 
T. Honma, Y. Haga, E. Yamamoto, Y. Onuki, H. Yam- 
agami, S. Takashi, and T. Yanagisawa, Phil. Mag. 79, 
1045 (1999). 

® D. Aoki, G. Knebel, I. Sheikin, E. Hassinger, L. Malone, 
T. D. Matsuda, and J. Flouquet, J. Phys. Soc. Jpn. 81, 
074715 (2012). 

® R. Okazaki, T. Shibauchi, H. J. Shi, Y. Haga, T. Matsuda, 

E. Yamamoto, Y. Onuki, H. Ikeda, and Y. Matusuda, 
Science 331, 439 (2011). 

S. Tonegawa, K. Hashimoto, K. Ikada, Y.-H. Lin, 

H. Shishido, Y. Haga, T. D. Matsuda, E. Yamamoto, 
Y. Onuki, H. Ikeda, Y. Matsuda, and T. Shibauchi, Phys. 
Rev. Lett. 109, 036401 (2012). 

S. Tonegawa, S. Kasahara, T. Fukuda, K. Sugimoto, N. Ya- 
suda, Y. Tsuruhara, D. Watanabe, Y. Mizukami, Y. Haga, 

T. D. Matsuda, E. Yamamoto, Y. Onuki, H. Ikeda, Y. Mat¬ 
suda, and T. Shibauchi, Nature Communications 5, 4188 
(2014). 

T. Shibauchi, H. Ikeda, and Y. Matsuda, Phil. Mag. 94, 
3747 (2014). 

S. Takagi, S. Ishihara, M. Yokoyama, and H. Amitsuka, J. 
Phys. Soc. Jpn. 81, 114710 (2012). 

I. Kawasaki, 1. Watanabe, A. Hillier, and D. Aoki, J. Phys. 
Soc. Jpn. 83, 094720 (2014). 

P. Thalmeier and T. Takimoto, Phys. Rev. B 83, 165110 

( 2011 ). 

P. Thalmeier, T. Takimoto, and H. Ikeda, Phil. Mag. 32- 
33, 3863 (2014). 

H. Ikeda, M.-T. Suzuki, R. Arita, T. Takimoto, 

T. Shibauchi, and Y. Matsuda, Nature Physics 8, 528 

( 2012 ). 

J. G. Rau and H.-Y. Kee, Phys. Rev. B 85, 245112 (2012). 
A. Akbari and P. Thalmeier, Phys. Rev. B 90, 224511 
(2014). 

W. K. Park, P. H. Tobash, F. Running, E. D. Bauer, J. L. 
Sarrao, J. D. Thompson, and L. H. Greene, Phys. Rev. 
Lett. 108, 246403 (2012). 

F. Bourdarot, E. Hassinger, S. Raymond, D. Aoki, V. Tau- 
four, L.-P. Regnault, and J. Flouquet, J. Phys. Soc. Jpn. 
79, 064719 (2010). 

F. Bourdarot, S. Raymond, and L.-P. Regnault, Phil. Mag. 
94, 3702 (2014). 

C. Stock, C. Broholm, Y. Zhao, F. Demmel, H. J. Kang, 

K. C. Rule, and C. Petrovic, Phys. Rev. Lett. 109, 167207 

( 2012 ). 

O. Stockert, J. Arndt, E. Faulhaber, C. Geibel, H. S. 
Jeevan, S. Kirchner, M. Loewenhaupt, K. Schmalzl, 
W. Schmidt, Q. Si, and F. Steglich, Nature Physics 7, 
119 (2011). 

I. Eremin, G. Zwicknagl, P. Thalmeier, and P. Fulde, Phys. 
Rev. Lett. 101, 187001 (2008). 

A. Akbari and P. Thalmeier, Phys. Rev. B 86, 134516 

( 2012 ). 

S. Chi, D. T. Adroja, T. Guidi, R. Bewley, S. Li, J. Zhao, 

J. W. Lynn, C. M. Brown, Y. Qiu, G. F. Chen, J. L. Lou, 
N. L. Wang, and P. Dai, Phys. Rev. Lett. 101, 217002 
(2008). 

A. Akbari, 1. Eremin, P. Thalmeier, and P. Fulde, Phys. 
Rev. B 80, 100504(R) (2009). 

K. S. Nemkovski, J.-M. Mignot, P. A. Alekseev, A. S. 
Ivanov, E. V. Nefeodova, A. V. Rybina, L.-P. Regnault, 


F. Iga, and T. Takabatake, Phys. Rev. Lett. 99, 137204 
(2007). 

G. Friemel, Y. Li, A. Dukhnenko, N. Shitsevalova, 
N. Sluchanko, A. Ivanov, V. Filipov, B. Keimer, and D. In- 
osov. Nature Gommunications 3, 830 (2012). 

A. Akbari, P. Thalmeier, and P. Fulde, Phys. Rev. Lett. 
102, 106402 (2009). 

A. Akbari and P. Thalmeier, Phys. Rev. Lett. 108, 146403 

( 2012 ). 

P. S. Riseborough, Phys. Rev. B 45, 13984 (1992). 

P. Thalmeier and A. Akbari, “Proceedings of the 50th 
karpacz winter school: Quantum criticality in condensed 
matter,” (World Scientihc, Singapore, 2015) p. 44 
(arXiv:1412.2523). 

C. Broholm, J. K. Kjems, W. J. L. Buyers, P. Matthews, 
T. T. M. Palstra, A. A. Menovsky, and J. A. Mydosh, 
Phys. Rev. Lett. 58, 1467 (1987). 

C. Broholm, H. Lin, P. T. Matthews, T. E. Mason, W. J. L. 
Buyers, M. F. Collins, A. A. Menovsky, J. A. Mydosh, and 
J. K. Kjems, Phys. Rev. B 43, 12809 (1991). 

A. V. Balatsky, A. Chantis, H. P. Dahal, D. Parker, and 

J. X.Zhu, Phys. Rev. B 79, 214413 (2009). 

A. Hewson, The Kondo problem to heavy fermions (Cam¬ 
bridge University Press, 1993). 

P. Aynajian, E. H. da Silva Neto, C. V. Parker, Y. Huang, 
A. Pasupathy, J. Mydosh, and A. Yazdani, Proc. Nat. 
Acad. Sci. 107, 10383 (2010). 

T. Takimoto and P. Thalmeier, Phys. Rev. B 77, 045105 
(2008). 

‘‘1 M. B. Maple, J. W. Chen, Y. Dalichaouch, T. Kohara, 
C. Rossel, M. S. Torikachvili, M. W. McElfresh, and J. D. 
Thompson, Phys. Rev. Lett. 56, 185 (1986). 

A. Villaume, F. Bourdarot, E. Hassinger, S. Raymond, 
V. Taufour, D. Aoki, and J. Flouquet, Phys. Rev. B 78, 
012504 (2008). 

N. P. Butch, M. E. Manley, J. R. Jeffries, M. Janoschek, 

K. Huang, M. B. Maple, A. H. Said, B. M. Leu, and J. W. 
Lynn, Phys. Rev. B 91, 035128 (2015). 

T. T. M. Palstra, A. A. Menovsky, J. van den Berg, A. J. 
Dirkmaat, P. H. Kes, G. J. Nieuwenhuys, and J. A. My¬ 
dosh, Phys. Rev. Lett. 55, 2727 (1985). 

Y. Kohori, K. Matsuda, and T. Kohara, J. Phys. Soc. Jpn. 
65, 1083 (1996). 

N. Emi, R. Hamabata, D.Nakayama, T. Miki, T. Koyama, 
K. Ueda, T. Mito, Y. Kohori, Y. Matsumoto, Y. Haga, 
E. Yamamoto, Z. Fisk, and N. Tsujii, J. Phys. Soc. Jpn. 
(2015). 

Y. Kasahara, H. Shishido, T. Shibauchi, Y. Haga, T. D. 
Matsuda, Y. Onuki, and Y. Matsuda, New Journal of 
Physics 11, 055061 (2009). 

K. Yano, T. Sakakibara, T. Tayama, M. Yokoyama, 

H. Amitsuka, Y. Homma, P. Miranovic, M. Ichioka, 
Y. Tsutsumi, and K. Machida, Phys. Rev. Lett. 100, 
017004 (2008). 

E. R. Schemm, R. E. Baumbach, P. H. Tobash, F. Ron- 
ning, E. D. Bauer, and A. Kapitulnik, Phys. Rev. B 91, 
140506(R) (2015). 

E. Hassinger, D. Aoki, F. Bourdarot, G. Knebel, V. Tau¬ 
four, S. Raymond, A. Villaume, and J. Flouquet, 
arXiv:0909.4188 (2009). 



